The Unusual Superconducting State of Underdoped Cuprates 
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There is increasing experimental evidence that the superconducting energy gap Ao in the under- 
doped cuprates is independent ol doping concentration x while the superfluid density is linear in 
X. We show that under these conditions, thermal excitation of the quasiparticles is very effective in 
destroying the superconducting state, so that Tc is proportional to a; Ao and part of the gap structure 
remains in the normal state. We then estimate Hc2 and predict it to be proportional to x'^ . We also 
discuss to what extent the assumptions that go into the quasiparticle description can be derived in 
the U(l) and SU(2) formulations of the t-J model. 

PACS numbers: 74.25.Jb,79.60.-i,71.27.+a 

While the anomalous properties of cuprate supercon- 
ductors have been discussed from the very beginning, 
it has become clear in the past several years that it 
is in the underdoped region that the cuprates deviate 
most strongly from conventional materials, both in the 
normal and superconducting states. NMR,^ neutron^ 
and c-axis optical conductivity'^ indicate the existence 
of a pseudogap in the normal state, and photoemission 
experiments^'^ reveal that the pseudogap is of the same 
size and k dependence as the superconducting gap. Fur- 
thermore, as the doping concentrating x is reduced from 
optimal doping, the superconducting gap is constant or 
may be slightly increasing, while the transition temper- 
ature Tc is reduced. Thus a strong deviation from the 
BCS ratio between energy gap and Tc is to be expected. 
At the same time, London penetration depth and op- 
tical conductivity^ show that the Drude weight in the 
normal state as well as the superfluid density in the su- 
perconducting state are proportional to x, and a linear 
relation between Tc and the superfluid density have been 
noted. ^ The small superfluid density has led a number 
of authors to suggest that phase fluctuation may be the 
determining factor of Tc and that the pseudogap may be 
interpreted as superconducting fluctuations.*'^ A related 
interpretation in terms of Bose condensation of pairs has 
also been suggested. ^'-^° A second school of thought starts 
with strongly correlated models such as the t-J model 
and interprets the pseudogap as the spin excitation gap 
in some RVB singlet state. In particular, this scenario is 
realized in a decomposition of the electron into a fermion 
which carries spin and a boson which represents a va- 
cancy in order to enforce the constraint of no double 
occupation of the t-J model. At the mean field level, 
spin and charge are separated, and in the underdoped 
region of the mean field phase diagram the fermions are 
paired in some intermediate temperature range, and be- 
come a d-wave superconductor only below the Bose con- 
densation temperature of the bosons. ^^'^^ Very recently, 
a modification of this mean field theory was proposed. 



which incorporates an SU(2) symmetry known to be im- 
portant at half filling. -'^^ It was argued that the new SU(2) 
formulation allows a smoother connection to half-filling 
and includes low lying fluctuations ignored in the origi- 
nal U(l) formulation. Features in the photoemission data 
are qualitatively explained by this approach, but so far 
discussions have been limited to the normal state. 

In this paper we examine the superconducting state of 
the underdoped cuprates. We begin with a phenomeno- 
logical approach, based on the existence of well defined 
quasiparticles in the superconducting states. We show 
that the combination of a large energy gap and small 
superfiuid density leads to unusual features. A key re- 
sult is that the superconducting state is destroyed by the 
thermal excitation of just the low lying quasiparticles, 
leaving the large energy gap intact. We believe this is a 
more effective mechanism of destroying superconductiv- 
ity than the phase fluctuation scenario. With very few 
assumptions we derive expressions for the temperature 
dependence of ps and for Tc and Hc2- We then consider 
whether the quasiparticle approach can be derived from 
the microscopic treatments of the t-J model. We con- 
clude that the U(l) formulation fails to obtain the correct 
temperature dependence of ps-, even if gauge fluctuations 
are included in the Gaussian approximation. We then in- 
dicate how the assumptions of the phenomenology can be 
derived from the SU(2) formulation, provided that quan- 
tum fluctuations of low lying excitations are taken into 
account. 

We limit our discussion to clean superconductors and 
assume that the elementary excitations in the super- 
conducting state are well defined quasiparticles with 

dispersion E{k) — f {sj^ — pY + ^1 ) where A^ — 

^Ao(cosfca;a -f coskyo) is a d-wave gap with a maxi- 
mum of Aq at (0,7r). In a tight binding parametriza- 
tion, e^ ~ 2tf{coskxa -\- coskyu). There are 4 nodal 
points. In the vicinity of the nodes near (^, ^j, we have 
the anisotropic Dirac spectrum E{k) = {v'jkf -\- w|A:|)^/^ 



where fci — {k^ + ky — 7r/a) /\/2, ^2 = (^x ~ ky)/V2, 
V2 = -y^Aoa, and Vf « 2V2/:/a for /i w 0. We now make 
the key assumption that in the presence of an electromag- 
netic field the quasiparticlc spectrum is shifted according 
to 



E{k,A) ^E{k) 



:^k-^ 



(1) 



where vi, = dsk/dk is the normal state velocity. Equa- 
tion (1) is correct to first order in A in the BCS theory. 
It has the consequence that the current carried by the 
quasiparticlc is —cdE/dA = —evk, i.e. it is the same as 
the current in the normal state and is different from the 
group velocity dEk/dk. The difference arises from the 
fact that the superconducting quasiparticlc is a superpo- 
sition of particle and hole states. For high Tc supercon- 
ductivity, we do not have the equivalence of BCS theory, 
so Eq. (1) must be regarded as an assumption, albeit 
a very reasonable one. This is particularly so near the 
node, where it is reasonable to believe that at the node 
the superconducting quasiparticlc should be the same as 
the normal quasiparticlc, so that its current should be 
given by t;;^. 

We next argue that the superfluid tensor defined by 
• _ efV^ gg^]^ i^g written as — pf,„ — —p^iT = 
0)(5^i/ — ^P^v where p'^{T = 0)/m = x/a^m is directly 
measured by the weight of the Drude peak in the nor- 
mal state and by A^ where Xl is the Landau pene- 
tration depth in the superconducting state. By tak- 
ing Xl ^ I6OOA for YBCO6.95 and x = 0.2, we find 
m = 2.lTOe. It is convenient to fit this mass to the bot- 
tom of a tight binding bandwidth with hopping integral 
tfi, and we have tf^ — {2ma^)^^ = 0.122 eV which hap- 
pens to be very close to J. On the other hand, -^p?ti, is 
given by the quasiparticlc response to the A field. This 
we can calculate by writing the free energy in terms of 
non-interacting quasiparticles, i.e., 



FiA, T) = -kTY^^n (1 + e-^^(^^^ 



fc.o 



(2) 



and differentiating twice with respect to A. We note that 
the neglect of quasiparticlc interaction is justified in the 
limit of small T and A because the density of states of 
quasiparticles vanishes linearly with energy, in contrast 
to the case of a Fermi liquid. As explained by Leggctt,^'* 
the Landau parameters enter in the form of mean field 
theory in Fermi liquid theory. The vanishing of the free 
quasiparticlc response functions implies that the Landau 
parameters play no role in this limit. We therefore obtain 



1 



/J.L' 



— 2E 



dE dE df 



dA„ dA„ dE 



(3) 



Strictly speaking, there is an additional term of the form 
2Y,ki(^'^E/dA^,dAy)f{E). This term vanishes in BCS 



theory due to particle-hole symmetry and we shall as- 
sume that it is also negligible in the present case. Using 
Eq. (1), we replace dE/dA^j, by the normal state velocity 
fp and ps (T) can be evaluated in a straightforward way 



m 



— aT 



(4) 



where a = {2en{2)/n)vF/v2 = (8fe(2)/7r)t//Ao. We see 
that for small x, the quasiparticlc excitation is an ef- 
fective way of destroying the superconducting state by 
driving p* to zero. By extrapolating Eq. (3) to p" = 0, 
we can estimate Tr as 



kT,^0.88xAo{th/tf) 



(5) 



If we assume that Aq is independent of x for underdoped 
cuprates, we see that Tc is proportional to x (or more 
precisely to Ps{T — 0)/m), thus providing an explana- 
tion of Uemura's plot.^ We shall see that tf^/tf w 1.8, so 
that the absolute value of T^ given by Eq. (4) is not un- 
reasonable. In particular, it is lower than the estimates 
based on phase fluctuation or Bose condensation, which 
typically gives Tc ~ xth- We emphasize that our mecha- 
nism is completely different from these other pictures, in 
that the quasiparticlc spectrum and the energy gap Aq 
comes into play. Obviously, Eq. (4) implies a strong de- 
viation from the BCS ratio between Tc and Aq for small 

X. 

Another important implication is that superconductiv- 
ity is destroyed when only a small fraction of the quasi- 
particles (with energy < xAq) are thermally excited. 
Thus the gap near (0, tt) must remain intact in the normal 
state, leaving a strip of thermal excitations which extend 
a distance proportional to x from the nodal points. This 
is qualitatively in agreement with the photocmission ex- 
periment. Of course our phcnomcnological picture does 
not provide a description of the normal state. It simply 
states that the normal state gap is an inescapable conse- 
quence of a finite Aq and a vanishingly small superfluid 
density as a; ^- 0. 

The fact that dps/dT is independent of a; and that both 
Ps and Tc arc proportional to x means that a scaled plot 
of ps{T)/ps{0) vs T/Tc should be independent of x for 
small T/Tc- In fact, such a scaled plot for YBCOg.gs and 
YBCOg.eo shows a remarkable universality over the entire 
temperature range. ^^ We can use the data to extract the 
ratio vp /v2 or tf/Ao using Eq. (3). Using the YBCO6.95 
data, we obtain a velocity anisotropy wf/w2 = 6.8. A 
slightly smaller ratio (by about 15%) as obtained from 
the YBCOe.eo data. This implies that tf/Ao = 1-7. If 
we assume Aq — 40 meV, we find tf = 68 meV, so that 
th/tf = 1.8 as mentioned earlier. Our value of tf im- 
plies a half-filled bandwidth ^ Atf w 270 mcV which 
is consistent with the photocmission data. This gives 
Vf = 1.18 X 10^ cm/sec. 



Equation (1) implies that in the presence of a magnetic 
field, the quasiparticle spectrum is shifted so that some 
of the quasiparticles are occupied in the ground state and 
a finite density of states is generated at the Fermi energy. 
This effect was predicted by Volovik^^ and observed in 
specific heat measurements.^^ We shall now use a similar 
picture to estimate the size of the vortex core and esti- 
mate Hc2. The idea is that as we approach the vortex 
core, according to Eq. (2), the local superfluid density is 
reduced. We identify the core size as the point when ps is 
reduced to zero, i.e., when the critical current is reached. 
Let us approach the core in the x direction. The A field 
is in the y direction and 2eA/ c should be replaced by the 
gauge invariant current V6' — —A. Near the vortex core 
the V6' term dominate so that we can approximate eA/c 
by (2i?)-iy. At T = we combine Eq. (1) and (2) and 
use —df/dE = S{E) to show that the local normal fluid 
density at the point Rx is given by 



^{Rx) 
m 






(6) 



where N{E) — ^^ 6 {E — E{k)) is the quasiparticle den- 
sity of states for one node. It is easy to show that 
N{E) = El{2-KV}V2)- By setting Eq. (5) equal to 
/9''(r = 0)/?Ti, we estimate the core size in the x direction 
to be 



Ri 



1 



(w/Z^Ao) [tfl^th 



(7) 



Note that it is greater than the BCS coherence length 
w/ZttAq by the factor x^^ . On the other hand, using the 
Tc estimate in Eq. (4), we can write 

Ri ^ VF/{l.6nkTc) . (8) 

which is quite close to the BCS coherence length written 
in terms of Tc instead of Aq. The two ways of writing 
the coherence length are of course equivalent in BCS the- 
ory, but very different for underdoped cuprates. A main 
conclusion of this paper is that Eq. (8) is the proper 
expression for the coherence length. 

When the core is approached from the (1,1) direction, 
similar considerations show that the core size is given by 
\/2Ri. Thus the core takes on an approximately square 
shape. This anisotropy in fact extends over long dis- 
tances and is responsible for the tendency to form square 
vortices in d-wave superconductors.^* We estimate Hc2 
by assuming that the square vortex cores are closed- 
packed, so that 



Hc2 - {hc/2e)/ml 



(9) 



Using Eq. (7) and the parameters we extract from the 
penetration depth data and setting Aq — 40 meV, x 
= 0.125 for YBCOe.e, we estimate i?i « 18.5 A and 
Hc2 ~ 146 T. Due to the crudeness of the extrapola- 
tion process, we expect both the Tc and Hc2 expressions 



to be overestimates. Perhaps a more realistic estimate 
may be made by using Eq. (8) for i?i and using the 
experimental Tc and we obtain instead i?i f* 30 A and 
Hc2 ~ 56 T which is close to the measured value of 40T.^^ 
While the absolute value of Hc2 is quite uncertain, the 
prediction that Hc2 is proportional to x^ (or more ac- 
curately to pliT = 0)), as long as Ao is constant for 
underdoped cuprates, should be amenable to experimen- 
tal test. The ideal systems to test this correlation are 
underdoped YBCO or Hg cuprates, which fall on the Ue- 
mura plot so that PsiT = 0) is proportional to Tc and can 
be accurately determined. In principle LACO is a good 
testing ground because x can be varied. Unfortunately, 
there are serious disorder effects for x < 0.1 and the nom- 
inally pure compound x ~ 0.15 is off the Uemura plot, 
for reasons which are not presently understood. Equa- 
tion (1) breaks down in the presence of disorder, and the 
low lying excitations may even be localized,^" so that the 
results in this paper are restricted to the clean limit. 

Next we comment on whether existing microscopic 
models can reproduce the assumptions of the quasiparti- 
cle description. In the U(l) formulation of the t-J model 
the normal state in the underdoped limit is described by 
d-wave pairing of fermions, so that there exists an en- 
ergy gap Aq which remains finite as x ^ in the normal 
state. ^^'^^ Superconductivity is driven by condensation of 
bosons and well defined quasiparticles are formed. The 
superconducting Tc occurs as an energy scale of A^nxth at 
the mean field level, and may be suppressed by gauge 
fluctuations.^^ In this theory the superfluid density is 
given by the loffe-Larkin rule,^^ ^^(7") ^ p^^ + pl^^ . 
Since the energy gap appears in the fermion spectrum, 
we expect p%{T) = {\ — x) — T/Aq while p^ k, x with 
a higher power in T correction. Then the U(l) theory 
predicts p'{T) = x- x'^T/Aq. While the T = value is 
correctly given to be x, the temperature dependence is in 
strong disagreement with Eq. (3) and with experiment 
in that a is suppressed by x^ . The origin of this difficulty 
is that the fermion does not couple directly to A, but to 
the U(l) gauge field a while the bosons couple to A + a. 
The external A produces a finite a but its magnitude is 
reduced by x. In the quasiparticle language, the shift of 
the spectrum in the presence of A is smaller than that 
given in Eq. (1) by x. It is difficult to escape from this 
conclusion in the U(l) theory, because gauge fluctuations 
are included at the Gaussian level in the loffe-Larkin rule, 
which should be a good approximation in the supercon- 
ducting state in that the gauge mode is strongly gapped 
by the Anderson-Higgs mechanism. 

It was shown recently^'' that the U(l) formulation does 
not connect smoothly to the half-filled limit which is 
known to exhibit an SU(2) symmetry. For small x, there 
are indeed low lying gauge fiuctuations with energy scale 
of order xAo which arc ignored in the U(l) formulation. 
A new SU(2) formulation was introduced, which allows 



these low energy fluctuations to be described in a natu- 
ral way. The low energy effective theory contains a boson 



part and a 
see Ref. 13 


fermion part: ie// = Lt + Lf. (For 
) The boson part is given by 
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Lb = 


- ibHdt 


— ieAo — iaoT 
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{blby-^iblb- 


-D2 


^{\b^\' ^ \b,\r 


(10) 



where b 



bi 



, and Di_2 are order one coefficients. 



The fermions arc in a staggered flux phase and couple 
only to the a^ gauge field. When |6i| = I62I 7^ 0, the sys- 
tem is in a superconducting state (which corresponds to 
the d-wave paired state in the U(l) formulation). When 
61 7^ and 62 = 0, the system is in a metallic state 
(which corresponds to the staggered flux phase in the 
U(l) formulation). Since D2 > the ground state is 
the superconducting state. The normal state at finite 
temperatures contains no boson condensation and is a 
state which fluctuates between d-wave pairing and the 
staggered flux phase of fermions. The fermion spectrum 
acquires a gap Ag which is finite for small x. 



k^Bosory', 
i/ \-^ \a 



terparticle spacing x~^/^. Secondly, Hc2 ~ x^ is a weak 
field in the sense that when compared with the hole den- 
sity x, the number of Landau levels occupied is x^^ 3> 1 
so that we are outside of the quantum Hall limit. Thus 
we expect the state for H > Hc2 to be a metallic state 
and the key question is what kind of metallic state it is. 
It is clear from the present discussion that for x <C 1, the 
energy gap at (0, tt) survives inside the vortex core and 
therefore for H > Hc2- The magnetic field drives a region 
of gapless excitations in the Brillouin zone which extends 
a distance x from the nodal positions, qualitatively sim- 
ilar to the normal state above Tc- It seems to us that 
two possibilities remain. First the gapless excitations are 
well defined quasiparticles in the Landau sense. In this 
case, Luttinger theorem requires that a breaking of trans- 
lation symmetry must occur to produce the energy gap 
and the metallic state may be understood as some form 
of staggered flux phase. Alternatively, the gapless excita- 
tions are not Landau quasiparticles, but acquire residual 
width due to quantum fluctuations, making this state a 
genuine non-Fermi liquid state. This latter scenario is an 
exciting possibility which deserves further investigation. 
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FIG. 1. Coupling between Ai and fermion quasiparticles. 

We note that, in the superconducting state (|6i| — 
1 62 1), ttfj, and Afj, decouple under the meanfield approxi- 
mation. Therefore the low lying fermion quasiparticles do 
not couple to the external electromagnetic gauge field A^ , 
and cannot reduce the superfluid density within mean- 
field theory. However unlike the U(l) case, quantum fluc- 
tuations of the gauge fields and the quantum fluctuations 
between the two bosons (both are omitted in the U(l) 
formulation) are important even at T = 0. Those quan- 
tum fluctuations induce a coupling between the fermion 
quasiparticles and the gauge potential A. One such con- 
tributions is illustrated in Fig. 1. We find that the shift 
of the quasiparticle spectrum is of the form given by Eq. 
(1), except that the Vk ■ A is multiplied by a numerical 
constant of order unity. Thus the SU(2) theory incorpo- 
rates the main ingredients underlying the present paper, 
i.e., a finite gap Aq, a superfluid density proportional to 
X, and a quasiparticle spectrum given by Eq. (1). Details 
of this microscopic theory will be given elsewhere. 

We believe that our prediction that Hc2 ~ x^ is sig- 
nificant for two reasons. First, it is in contrast with 
models based on Bose condensation which should predict 
Hc2 ~ X since the coherence length in that case is the in- 
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